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We investigate the evolution of a population of non-interacting particles which undergo diffusion 
and multiplication. Diffusion is assumed to be homogeneous, while multiplication proceeds with 
different rates reflecting the distribution of nutrients. We focus on the situation where the distribu- 
tion of nutrients is a stationary quenched random variable, and show that the population exhibits 
a super-exponential growth whenever the nutrient distribution is unbounded. We elucidate a huge 
difference between the average and typical asymptotic growths and emphasize the role played by the 
spatial correlations in the nutrient distribution. 
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I. INTRODUCTION 



The evolution of a population in an inhomogeneous environment with spatially varying growth rates can display 
complex dynamical patterns, resulting from the competition between diffusion, random multiplication and possibly 
advection [SIl]. Inhomogeneities of the environment can greatly affect the population dynamics resulting in anomalous 
spreading (neither diffusive, nor ballistic) and intermittent behaviors (for a review on intermittency in random media 
see, e.g., [3]) with patches of population growing in favorable places and being surrounded by large desertic regions. 
Diffusion tails play the role of vanguard agents that explore hostile regions, in search of ever-more hospitable locations 
to settle and to multiply. The offspring in the new settlement will eventually outgrow the original colony and generate 
a strong gradient current: this can be seen as a migration of the whole population. A simple and concrete example 
of this problem is provided by bacteria multiplying and diffusing in a petri dish where nutrients and inhibitors are 
unevenly distributed, resulting in complex growth patterns that are observable in actual experiments [11 16] . In a 
more abstract setting, the actual inhomogeneous space can be replaced by a rough phenotypic landscape in which the 
fitness functional (that governs the growth rate) takes different values: here the coordinate variable does not label a 
spatial location but rather the genetic content of each individual. This point of view has been inspired by S. Wright, 
R. A. Fisher and J. B. S. Haldane's classical studies in population genetic [5]. More recently it was developed by M. 
Eigen [7] as a model for punctuated evolution of quasi-species and has been further investigated by W. Ebeling et al. 
[5], Y.-C. Zhang [3], M. N. Rosenbluth [T^ and many other authors; see [UlIIl] for a review of more recent work. 

In the present work, we examine a population of non-interacting particles that diffuse and undergo the birth/death 
process depending on the availability of nutrients. The population density n(x, t) evolves according to the diffusion 
equation with a multiplicative noise: 

9nfx t) 

— ^ D\/'^n(x,t) +r](x)n(x,t) with n(x, 0) no(x). (1) 
dt 

The diffusion coefficient D is assumed to be uniform, while the birth/death rate 7/(x) is inhomogeneous. We shall 
focus on the simplest case of the stationary noise, rj = ?7(x), and investigate how spatial correlations of the noise affect 
the growth rate of n(x, t). 

The Langevin equation ([l]) is linear, reflecting the basic assumption that particles do not interact. The noise ?7(x) is 
known through its stochastic properties and changing these properties may drastically affect the behavior of n(x, t). 
Equation ([I]) can be generalized in various ways, e.g. one can take into account advection |13H16j . study the evolution 
of a vector field in a random background (e.g. a magnetic field in the dynamo effect) [2], or consider several coupled 
fields as in the case of chemotaxis [T71 |lH| ■ One can introduce a non-linearity in order to take into account saturation 
effects [in]. Besides, it can also be interesting to consider situations with rj = 7/(x, t), the dependence on time reflecting 
e.g. seasonal variations. 

Although we shall chiefly employ a population dynamics vocabulary (particles, nutrients, migrations), it is useful to 
keep in mind that the stochastic partial differential equation ([T]) arises in many different contexts. In chemical physics, 
reaction kinetics is often modeled by equations similar to Eq. (ij and these equations allow to predict macroscopic 
patterns in the spatial distribution of reagents [5DJ[2T]. When the amplification rate r](x) takes only negative values, 
equation ([I]) describes non-interacting particles that diffuse in a medium with random absorption. In the extreme 
case of rj{x) = — '}2n 5{x — Xi) with random positions Xi of the traps, the trapping sites become perfect because n[x, t) 
has to vanish for x = x^. (More precisely, the above formulation applies in one dimension; in higher dimensions, 
the traps should be finite and absorbing conditions are set on the boundaries of the traps.) The literature on this 
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classical subject is vast; one of the most celebrated results asserts that the density decreases according to a stretched 
exponential law: In n ~ — f^/ ('^+^) in d dimensions J22f[29] . 

Another important interpretation of equation (IT]) in terms of polymer dynamics can be obtained from its formal 
solution, 



7i(x, t) = j dyK{^, t; y, 0) no(y) , (2) 
where the Green function K{x,t]y,0) is given by the Feynman-Kac formula |30| 



is:(x,t;y,0) = Ey^x exp 



?7(x(r))(iT 



(3) 



The expectation value Ey_>.x with respect to the Wiener measure is taken over all paths x(t) that begin at y at t = 
and end at x at time t. This formula can be rewritten in more familiar manner as a path integral 



(x, y, 0) = y ^ _ I?x(t) exp| dr ^- + r/(x(r))J | 



(4) 



It is then natural to interpret x(r) as a c?-dimensional Gaussian polymer in the random potential ?7(x) [3TI - I34] . or 
as a d + 1-dimcnsional directed polymer with columnar disorder [35 (for a review on directed polymers see |36|). 
The physics of directed polymers in a random medium is closely related to the growth of random surfaces. In one 
dimension, for instance, we write n — e"^ and recast ([I]) into rrit ~ D{m^ -'rrnxx) Differentiating this equation 

with respect to x and writing h — -~2Dmx we obtain 

ht + hhx = Dhxx + where ^ ~ —2Drjx ■ (5) 

Thus after the transformation the noise becomes additive, but the governing Eq. ^ is now a non-linear stochastic 
partial differential equation. The Langevin equation ^ formally resembles the Kardar-Parisi-Zhang (KPZ) equation 
[36l [37] , yet the noise ^ in equation (J5| is stationary, whereas in the KPZ equation the noise depends both on space 
and time. [In the KPZ equation, the noise is typically taken to be the Gaussian white noise, (?7kpz(x, t)?7KPz(x'i t')) = 
2r (5(x — x') (5(i — t').] Hence the physical properties of the solutions to Eq. (J5]) are significantly different from those 
of the standard KPZ equation [55] . 

Equation ([T]) can also be viewed as a Schrodinger equation in imaginary time and with a random potential; it 
is thus related to the physics of locahzation [521 13H]- In the present study n{x,t) is a density of particles (not a 
wave-function). For instance, in the no-noise case the total mass J dxn{x,t) is conserved, whereas the analogous 
integral in the quantum case is not conserved. Furthermore, we are not directly interested in eigenstates of (|T]) but 
rather in the temporal behavior of its solution starting from a localized initial state. All these differences explain why 
the phenomenology is rather different from that of localization, although some of the techniques developed to study 
quantum disordered systems are useful in the analysis of Eq. ([T]) . 

The behavior of solutions of the apparently simple linear Langevin equation ([T]) is not yet fully understood due to a 
number of puzzling features. One such feature is an astonishingly fast growth (and sometimes a blow up that occurs 
in a finite time, or even instantaneously |10j). There are a few causes of these striking behaviors: 

• The noise is multiplicative. 

• In many simple models, the noise is unbounded. Hence, the regions with large positive 77 play a dominant role 
and lead to a counter-intuitive super-exponential growth. In this situation, the discretized (in space) versions 
of Eq. ([1]) differ drastically from the strictly continuous version. 

• The lack of self-averaging which is manifested by the huge difference between average and typical behaviors. 

The goal of this work is to analyze the Langevin equation ([T]) when the noise is strongly correlated in space. Some 
of our results are presented in Table l] where we display only the scaling laws (numerical constants will be given 
throughout the text). The part of Table |l] describing the asymptotic growth laws for the uncorrelated Gaussian white 
noise summarizes previous work. The notion of "correlated" noise is of course a bit vague. The one-dimensional case 
is exceptional as the very natural assumption that the increments of the noise are uncorrelated leads to the correlated 
noise r]{x) which is essentially a trajectory of a random walk (in the lattice setting) or a Brownian motion (in a 
continuum setting), with x playing a role to time. Table |l] presents the asymptotic growth laws corresponding to such 
a noise. (Our main higher-dimensional results are collected in Table [h]) 
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Uncorrelated 


Correlated (d — 1) 


Lattice substrate 

e* (Average, all d) 
(Typical, all d) 


Random walk landscape 

e* (Average) 
e*^^'' (Typical) 


Continuum substrate 

e*' (d = 1) 
Blow up (d > 2) 


Brownian landscape 

e* (Average) 
e*'^^^ (Typical) 



TABLE I: Asymptotic growth laws for some of the population dynamics (numerical factors are omitted). The noise is assumed 
to be Gaussian, apart from the case when the noise forms a random walk landscape, that is, the increment of the noise are 
bounded. (The uncorrelated bounded noise is investigated in Sect. II B ) The results for the correlated noise in higher dimensions 
are discussed in Sect. IIII"C] 



The rest of this paper is organized as follows. In Sect. [TTj we review the behaviors when the noise is either 
uncorrelated or has short range correlations. We emphasize the emergence of puzzling behaviors and explain how 
seemingly contradictory results scattered in the literature can be synthesized in a coherent way. In Sect. |IIl| we analyze 
the much less studied situation of a correlated noise and determine the growth law for the total population size in the 
one-dimensional setting. We then qualitatively describe the situation in higher dimensions. 



II. POPULATION DYNAMICS WITH SHORT-RANGE CORRELATIONS 



In this section we assume that the local growth rate ?y(x) is on average homogeneous. Therefore, the mean 
value /i = (r/) is constant which can be set to zero; the general case is recovered by redefining the local density: 
n(x, t) — >■ n(x, t) e''*. Further, because of homogeneity, the spatial correlations of 77 must be translationally invariant. 
The simplest assumption, customary in studies of Langevin equations, is to consider the random potential rj to be 
uncorrelated at different spatial locations. The fluctuating potential 77(x) is thus taken to be a Gaussian white noise 
with zero average: 

(r/) = 0, (r7(x)7;(x')) ^ 2r5(x - x') . (6) 

The stochastic properties of the noise are now fully specified and equation ([I]) defines a well-posed problem which has 
been studied in numerous works [U |3l [8HI0I l34l [39'-'43' mostly in one dimension. The conclusions of these studies, 
established through various methods and approximations, seemed initially contradictory. However, the issue was 
settled in |10j by an exact analytical calculation based on a minimax variational principle. In the large time limit, the 
dominant contributions to the population density arise from very small, rapidly growing isolated regions. Although 
these regions are very rare and highly improbable, they produce high density peaks that dominate the whole statistics. 
This has led to the following unexpected asymptotic behavior in one dimension: 

\n{n{0,t))^^ . (7) 

This super-exponential growth has been also found by estimating Brownian motion expectations |40j and by a path- 
integral approach that has additionally allowed to calculate the pre- factors of the exponential behavior [i^ : 

p2.5/2 /T^t^\ 
W0,^))^;T7^T^exp — . (8) 



In higher dimensions, the behavior is even more puzzling. It has been argued 'IQ] U?] that in two dimensions, a 
divergence occurs at a finite time tc ~ D/T; for t > tc the solution blows up. Further, when d > 2, the divergence is 
instantaneous, i.e., equation ([T]) with Gaussian white noise (|6| is ill-defined. These striking behaviors exhibited by the 
Langevin equation ([T]) with Gaussian noise ^ are consequences of two major properties of the noise: The Gaussian 
noise is both uncorrelated and unbounded. We now investigate the consequences of relaxing these assumptions. 
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A. Taming the White Noise: Lattice Regularization or Finite Short-Range Correlations 



The white noise displays totally uncorrelated fluctuations at all scales. However, in the present context of a 
population dynamics it is natural to assume some spatial coherence in the local varying conditions and that the noise 
has a non-vanishing correlation length. The simplest manner to implement this model is to put the system on a lattice 
and to assume that the noise is uncorrelated at different lattice sites. The effective correlation length is thus equal to 
the lattice spacing. Mathematically, one has to solve the discrete set of equations 

nj = DV'ui + njT/j (9) 

where i = {ji, ■ ■ ■ ,jd) G '^'^ if the lattice is hyper-cubic. Further, the operator denotes the discrete Laplacian (e.g. 
in one dimension we have V^rij = '^-j-i ~ 27ij -I- The noise in ^ is Gaussian with the following characteristics: 

iv) = 0, iViVi) = 2r5i.j . (10) 

In this case, the growth is chiefly universal, it is independent of the spatial dimension (up to numerical constants) 
and is given by (n) ~ e^* [Tj. Hereinafter, we shall use the convention that in the asymptotic such as n ~ exp{At°'), 
the displayed term gives the correct controlling exponential factor, so that the actual asymptotic may be something 
like n ~ Bt^ exp{At°') . Thus the asymptotic (n) ^ e'"* is actually a shorthand formulation of the leading asymptotic 
of the logarithm: 

(11) 



The growth law (11) involves averaging over the disorder (12 1. The simplest set-up that automatically enforces 
such an averaging occurs in the situation when the evolution begins from the uniform initial condition: nj (0) = 1 for 
all j S Z'^. Indeed, no averaging is needed because, for the infinite lattice, all values of the noise are appropriately 
sampled. Thus, we can ignore diffusion altogether J]. Then n = e''*, so that 

eJ,t-'^]=e-^\ (12) 



/AnT V 4r 



If, however, the initial condition is localized, e.g. nj(0) = Syo, the behavior (11 1 will arise only after averaging over 
all distributions of the disorder. The typical behavior with a fixed noise, however, differs greatly: this is the sign of 
the lack of self-averaging. To establish the typical growth law, we denote by L the size of the domain visited by the 
particles. This domain contains about L'^ sites. The largest noise 77* (L) at these sites is evaluated using the extreme 
statistics criterion (see e.g. [IS]) to give 




from which 



(13) 



,(L) ~ y/4dT\nL. (14) 



We now ought to find out how the size L of the domain visited by a particle grows with time. A naive estimate 
postulates a diffusive scaling law L ~ y/i, thereby leading to rityp ^ exp{r]^{L)t) ~ exp(ty/2dr \nt) . This is wrong, 
however, as explained e.g. in Ref. [1]. The correct argument proceeds by averaging the optimal growth exp{ri^,(L)t) 
for a given path of lengths L > over all possible paths weighed by their probability of occurrence. Therefore, the 
typical population size grows as 

We calculate this integral by the saddle-point method. The exponent has a sharp maximum at 

^^(hJyVi' T^tD'/'iMTy/\ (16) 



Keeping only the dominant exponential factor in ( [15^ we obtain 

n,y,{t) ^ e'^^^ . (17) 
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This typical growth is essentially universal (the spatial dimensionality appears only in amplitudes) and the growth 
is just barely faster than exponential. It is important to note that the typical population growth rate results from 
particles that follow optimal paths that are almost ballistic, as seen from equation ( [T6| , rather than diffusive; these 
optimal paths are therefore highly non-typical individual trajectories. Here, the naive estimate L ^ \/t does, by chance, 
provide the correct functional dependence t\f\x\A, inside the exponential in equation (17), but with a coefficient wrong 



by a factor \f2. In the next section, we shall encounter cases where the naive estimate leads to erroneous results. 



The average population growth, given in equations (11) or (12), radically differs from the typical growth (17) 



However, these two results can be reconciled as follows. We have established (17) by estimating the value of rjJjT), 



the largest noise that occurs amongst L sites, using the criterion (|T3|. Yet the value given in (14) for r]^{L) is valid 



for a typical realization of the noise rj[x). In fact, the average growth of the population, given by (n), is dominated 
by highly non-typical realizations of the noise that must be taken into account: in order to calculate (n) correctly we 
have to let both the path and the background noise fluctuate. For L independent realizations of the Gaussian random 
variable rj, the cumulative distribution of the maximum M is given by 



Prob(77„,ax < M) 




,-'?'/4r 



(18) 



Taking the derivative of this expression, we find that the probability distribution of M 



Pl{M)^L 




-r,V4r 



(19) 



where the last expression, in which we have retained only the controlling exponential factor, is valid for large values 
of M. Since the average population at time t over a range of L sites grows as e^* , we obtain 



dL 



exp(-4lJt) 



/o {ATTDty/^ 
The asymptotic is evaluated using the saddle-point technique to yield 



dMe^*e-*^'/4r. 



(20) 



(21) 



in agreement with equation ( 12 ) which was obtained in the zero-dimensional case. We note that the diffusion constant 



D appears neither in the average behavior nor in the typical behavior (17): it affects only the sub-leading corrections. 

Thus the short-ranged correlations drastically modify the behavior of the solutions to Langevin equation ([l]) by 
regularizing the noise term. The short-range fluctuations that were responsible for the blow-up in dimensions d > 2 
are suppressed and, on the lattice, equation ([T]) is well defined in all dimensions. The growth of n{x, t) is the universal 
Gaussian law (11) that does not change with dimension. We finally note that another way to regularize equation ([T]) 
without discretizing space is to consider a colored Gaussian noise instead of a white noise. A frequently used example 
is the Gaussian Ornstein-Uhlenbeck noise with correlations given by {ri{x)r]{x')) = ^ exp(— |a; — x'|/^). This noise has 

exponentially decaying correlations with correlation length ^. For such a noise, the population grows as (n) ~ e 
(see [10]). This asymptotic is again independent of the diffusion constant and the dimensionality of space. 



B. Taming the White Noise: Bounded Noise Distributions 



The lack of the upper bound for the Gaussian noise is an obvious reason for the appearance of the faster-than- 
exponential growth found in ([t]), ([IT), and (17). For a bounded noise with ry < 77max, the growth cannot be faster than 



Interestingly, in most cases t 



le controlling factor is equal to e''"""' and the spatial dimensionality or details of 



the noise distribution (such as the behavior of the noise distribution function p{ri) in the proximity of ry < r^max) play 
a secondary role, namely they affect the pre-factor in the growth law. Let us look at this pre-factor. To appreciate 
its behavior, it suffices to analyze noise distributions with a finite number of different values of the noise. Without 
loss of generality we set the maximal noise to unity and write 



piv) pS{v - 1) + X! ~ 



(22) 
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with 

n 

P + = ai,...,a„<l. (23) 

i=l 

If the diffusion and the mult iphcat ion process occur on a one-dimensional lattice, the lattice can be thought to be an 
array of domains where the noise is maximal. Adjacent domains are separated by sites where the noise is smaller. 
The probability density of domains of length k is given by 

nfe = (i-p)V (24) 

where the factor accounts for k consecutive sites with maximal noise and the factor (1 — p)^ assures that the noise 



at the boundary sites is smaller than 1. (Using (24) one can compute the fraction of lattice with maximal noise to 
yield X]fe>i ^^fe = p as it should be.) 

Consider the simplest situation where the entire lattice is initially uniformly filled; ^^(0) = 1 for all j. The 
asymptotic behavior can be quantified by the average density 

1 ^ 

n(i)= hm - ^n,(t). (25) 

i=i 

We first observe that the average density has the trivial upper bound 

n < e* . (26) 

We now construct a lower bound for the average density. The idea is to consider the evolution on domains where 
the noise is maximal and use the absorbing boundary conditions on the ends of each domain. This is an obvious 
lower bound, yet we will see that it exhibits largely the same growth as the upper bound. To proceed, we make 
the assumption (to be confirmed a posteriori) that the chief asymptotic is actually provided by very long domains 
{k^ 1). For such domains we can replace the discrete diffusion equation by a continuous one and we need to solve 

^ = DV^n + n (27) 

on the interval < x < fc subject to the initial condition n{x, 0) = 1 and the absorbing boundary conditions 
n(0,t) = n{k,t) = for i > 0. In the long-time limit, the spatial distribution approaches the smallest eigenfunction 
of the Laplace-Dirichlet operator: n = f(t)sm{Trx/k). Plugging this into (27) we get n ~ gti^-D-R /fc ) gjj^z™^ Using 



Eqs. (|24|) and (|25|), we arrive at the estimate for the lower bound: 

/>oo 



k>l 



where we have replaced the sum by an integral because the asymptotic is dominated by the contribution of large 
domains. This integral can be calculated by the saddle-point method. One finds that the exponential term in the 

1/3 

integrand has a sharp maximum at fc* = [27r^Di/ ln(l/p)] . Keeping only the leading and sub-leading terms, we 
arrive at the lower bound 



> expjt- ^ (27r2i:)i)i/3[in(l/p)]i/3| 



(29) 



Comparing the upper and lower bounds, Eqs. (26) and (291, we see that the controlling exponential factors are the 



same. This provides an evidence in favor of the general assertion that for an arbitrary bounded noise in arbitrary 
dimension the controlling exponential factor is universal and determined by the maximal noise: 

n ~ e""-''* . (30) 



The above derivation of the lower bound (291 can be generalized to an arbitrary dimension. We again consider 
domains of neighboring sites with maximal noise. In principle, there can be an infinite domain (when the density 
p of the maximal noise exceeds a percolation threshold Pc{d)). Let p < Pc{d): If in this situation the controlling 
exponential factor is still given by ( |30[ ), it will certainly be valid for larger p. When p < Pc{d), the domains are finite 
and generally small. We now proceed as before, namely we set n = outside the domains as this will obviously provide 
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a lower bound. A well-known argument (see |45j ) implies that the domains which lead to the largest contribution are 
balls, so one must solve 

inside the ball r < R with the absorbing boundary condition n(r = _R, <) = on its surface. The solution reads 
?^ ^ gtii->'iD/R ) ^ where Js is the Bessel function with index S — (2 ~ d)/2 and Ai is the first zero of this 

Bessel function. Proceeding as in one dimension one arrives at a lower bound 



n>e'^ dRR''-^expl^-Xl^-ln{l/p)VdR''^ (32) 

where we have taken into account the fact that the probability that all sites of a ball of radius R have maximal noise 
scales as p^'^^ for large R (here Vd is the volume of a unit ball). Estimating the integral we obtain 



n > e* exp{-Q(i.t)^/(-^)} with = '-±^ j^^^^^^MM) ' . (33) 

For the continuous noise distributions, the behavior near the maximal noise plays a crucial role. We have analyzed 
noise distributions that behave as 

Prob(l - e < ?7 < 1) = Ae", a > 1 (34) 



2/(d+2) 



in the e ^ limit. The corrections to the controlling exponential factor (30) are similar to the case of the discrete 
noise distributions, e.g. in one dimension 

n > e' exp{-C(a, A)ti/3(i^^)2/3| (35) 



III. POPULATION DYNAMICS IN A BROWNIAN LANDSCAPE 



In this section, we relax the unrealistic assumption that the noise is uncorrelated when the distance exceeds a 
certain threshold. Within the ecological interpretation where the noise refers to local conditions, it is natural to 
assume that conditions change from site to site, yet if somewhere conditions are very good, they are also very good in 
the proximity. This suggests to consider a model where ri{x) is a random landscape. In one dimension, one practical 
realization of such a landscape is to take r]{x) to be a random walk trajectory (in the lattice setting) or a Brownian 
trajectory (in a continuum setting). In higher dimensions, rj will be taken to be a Gaussian field. In all cases, the 
roughness of the surface defined by the noise governs the population growth law. 



A. One-Dimensional Case: Heuristic Analysis 

In one dimension, ri{x) is taken to be a Brownian curve, in which x plays the role of a 'time' variable. We assume 
that the initial population seed is located at the origin n{x, t — 0) = S{x) and we set ri{x = 0) = 0. Then the noise is 
given by 

r]{x) = / ^{u)du (36) 
Jo 

where ^ is a Gaussian white noise. Thus, the autocorrelation of the landscape reads 

Mxf) = 2r|x| . (37) 

In the following we normalize the Brownian landscape by replacing r]{x) by \/2rri{x). 

First, let us estimate the typical and the average growth laws by employing a heuristic reasoning. Heuristic 
arguments elucidate the physical mechanisms leading to the super-exponential growth and shed light on the crucial 
distinction between the average and typical behaviors. 



8 



To estimate the typical growth, one could argue that the particle visits roughly L ^ y/t different lattice sites (in 
one dimension) during the time interval (0, t) and that highest value of the Brownian noise among these sites is 



L ^ t^/^. Because the density must grow as e''"""=*, one could anticipate that n > 



. However, the above 



derivation is too rough even for a heuristic argument: as already discussed in the previous section, the Feynman- 
Kac formula Q, representing a formal solution to equation ([l]), contains a summation over all possible trajectories. 
Therefore we should let L fluctuate. Thus we write 



ntyp ~ exp 



tVL 



L2 
4t 



and maximize with respect to L to give L ^ t^^^ leading to 



,5/3 



Hyp 



(38) 



(39) 



Interestingly, the optimal length L ~ t^^^ is super-ballistic. 

Another argument leading to the same result proceeds by saying that, since the total number of particles grows 
very rapidly, the total number of visited sites actually grows faster than diffusively. If all particles which are present 
in the system at time t were initially at the origin (an admittedly rough assumption) we estimate the size L of the 



/4irt 



segment of sites visited by the particles from the criterion 

that n grows exponentially, i.e., much faster than a power law) L 
Vmax ^ V L ^ (t Inn)^/''. Now using n ^ g''"™* we obtain yy^ax 
'7max: wc get t (lu n) "'^^^ ^ Inn, leading to (39|. 



-LV(4i) 



' 1. This gives (we take into account 

The maximal noise on this segment is 
i^^ Inn. Combining these two expressions for 



The typical growth law (39) is obtained for a given realization of the potential r\{x). The average growth of the 



population, averaged over different realizations of the potential, is very different, namely it is much faster. Again, the 
very rare fluctuations of the landscape dominate the average. The simplest way to estimate the average growth is 
to keep two free parameters, the size L of the segment visited by the particles and the maximum M reached by the 
noise ^(x) on this segment. This leads to 



exp 



tM - 



2L 



1/ 
At 



(40) 



where the factor exp(— M^/2L) represents the tail of the distribution of the maximum of a Brownian path over a 
range L. The maximum distribution of a Brownian path is a classical result that can be derived via the image method 



Maximizing in L and M we get L t ,M ^ t and then Eq. (40 1 results in 



(41) 



The asymptotic growth laws ( 39 ) and ( 41 ) crucially depend on the assumption that the noise is the Brownian 



landscape. To illustrate other possible behaviors we give two examples. 



1. Random walk landscape 



In this situation, the population dynamics occurs on the one-dimensional lattice and rj is assumed to be a random 
walk rather than a Brownian motion: ?7j+i — rjj = ±1, where ± are chosen independently and with equal probabilities 
so that the surface ri{x) has no tilt. The noise can go up L steps in a row (albeit with a small probability 2~^) and 
this rare fluctuation provides the dominant contribution. Indeed, writing 



exp 



l2 

tL L In 2 

4< 



(42) 



we see that in the t — >■ oo limit the exponentially small probability of the rare fluctuation is totally outweighed by its 
huge contribution. Maximizing in L we obtain L ^ and then Eq. ( 42 ) leads to 



(43) 



The difference between (41 1 and (43) is the consequence of the fact that the noise increments are hounded for the 
random walk landscape. 
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2. Fractional Brownian motion 



Consider the case of a self-affine disordered landscape, in which Gaussian fluctuations grow with distance with a 
positive Hurst exponent < < 1, namely {'q{x)rj{y)) = \x — y\^^ ■ (The localization of a quantum particle in such 
self-affine potentials has been recently studied, see and references therein.) Then ?7inax ~ and the same 

reasoning as above leads to the growth law 

lnntyp~t(^+^)/(^-^). (44) 



B. One-Dimensional Case: Quantitative Analysis using the WKB Method 



In this subsection, we derive the asymptotics (39) and (41 1. Equation ([T]) is linear, so it is useful to perform a 
spectral decomposition. We write 



where the eigenfunction he^x) satisfies a Schrodinger equation 



EnE{x) = D 



dx'^ 



+ w2Tiri{x)nE{x) , 



(45) 



(46) 



with normalized Brownian landscape 'q[x) playing the role of a potential. We also impose the normalization condition 
on the total mass of the eigenfunction 



nE{x)dx 



(47) 



the eigenfunctions ns are mutually orthogonal, we obtain C{E) = ueIO) and this allows us to write 



The coefficient C{E) in equation (45 1 is determined by the initial condition. Using n{x,0) — S{x) and the fact that 

(48) 



n{x,t) = / dEe'^^nE{0)nE{x) 



Integrating over x and using the normalization ( 47 1 yields 



n{t) = J dEe'^^UEiO) . 



We now estimate n£;(0). First, we rewrite the Schrodinger equation (461 as 

d'^nE{x) 



dx-^ 



/2r 

~D~ 



{ri{x) ~£)nE{x) = 0, 



E 



E 



2T 



(49) 



(50) 



The Brownian landscape rj{x) starts at 77(2: = 0) = and it will remain smaller than £ up to a first crossing-point 
x*{£) > such that r/{x*) = £ (see Figure [l]). We know the typical scaling x*(£) ~ Because ri{x) is a non- 
monotonous function of x, there will be many subsequent crossing points. (For the Brownian landscape, there will be 
infinitely many crossing points immediately following the first crossing [49].) We are interested in the behavior in the 
vicinity of the origin, see Eq. (49), so these further crossing points play a little role. 

In the long-time limit, the dominant contribution into the integral in (49) is provided by large values of 'energy' E. 
(We shall confirm this assertion below.) Therefore we can analyze Eq. (50) using the WKB method [501151]. On the 
interval x < x* , the WKB solution reads 



nE{x) 



1 



{£~r,{x)f' 



exp 



-(2r) 



1/4 



D 



dx' - T](x') 



(51) 



We now define x' = E'^t and use the scaling property r\(EP'T) = £'I]{t). Then, from the expression (51), we obtain 
(keeping again only the dominant exponential factor) 



nsiO) exp - 



^5/2 
2TVD 



dr ^Br 



(52) 
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where r* represents the first moment when the Brownian motion = 1 — ri{T) that starts at 1 at r = 0, crosses the 
origin: B^- > for r < r* and Br* = 0. 

We now consider a typical reahzation rj of the landscape. The integral that appears in (52) has a constant value 

^5/2 



K,j that depends on the landscape. Thus 



"-£;,typ(0) exp -K, 



Substituting this expression in the spectral decomposition (|49|) leads to 

^5/2 



ntyp 



dE exp [Et-K 



exp 



3/4 



5 \bK,^ 



2/3 



(53) 



(54) 



The second asymptotic in (54 1 has been derived via the saddle-point method. The maximum of the integrand occurs 
at E' ~ i^/^ which diverges as i -> co] this justifies the use of the WKB solution (51 1. The final result (54| qualitatively 
agrees with the prediction of Eq. (39) which was established using heuristic arguments. 

We now determine the value of the average population {n) , where the average is taken over all possible realizations 
of the noise. First, we need to calculate 



exp 




(55) 



where the expectation value is taken over all Brownian paths starting at 1 at r = and vanishing at r* for the first 
time. This expression is the average of a first passage exponential functional of the Brownian motion and its value 
can be determined using the general method described e.g. in Ref. [52|. The procedure applies to a functional of the 
form 



Q{xo) = Q{xq;p,U) 



expi—p J dTU{Br) 



(56) 



that involves an arbitrary smooth function U{Bt-) replacing \/B^ which appears in (55). The normalized Brownian 
process starts at any a;o > and as a function of xq, the functional (56 ) satisfies the backward Fokker-Planck equation 



(57) 
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with the boundary conditions 

Q{xo = 0) = 1 and Q{xo +oo) = . (58) 



Our problem, equation ( 55 ) , perfectly fits into this framework. Thus we must solve the following differential equation: 



-——2 py^Q{xo) = where p= (59) 

z uXq zL \ U 



We need to calculate Q{1) because the Brownian paths in equation (55) start at xq = 1. The asymptotic behavior of 
the solution of Eq. ( [59| can again be obtained via the WKB method. Writing Q = e'^ we find S{xo) = — | ^y2p (xq)^^'^ 
in the leading order. Therefore 

/ 4^5/4 N 

(n^(0))~expU^^,^7^j . (60) 



Inserting this result into equation ( 49 1 and evaluating the integral by saddle-point method we arrive at 

(n) - exp{lDrh^) . (61) 



This exact asymptotic qualitatively agree with Eq. (41) which was derived using qualitative arguments. 

As a side remark, we notice that the growth laws in judiciously chosen deterministic landscapes ri{x) shed light on 
the growth laws in random landcapes. For instance, 

• If we consider 7j{x) — nx, then equation ([T]) admits an exact solution 

, , 1 / {x-nDt^f K^Dt^ 

n(x,tj = exp' ' 



yf^m \ ^Dt 3 

representing a Gaussian profile with uniformly accelerated center and with total mass growing as e*^. This 
solution can be obtained by the previously described spectral analysis: UEix) = ki[E — x) where Ai is the Airy 
function. This example can be used as a template for applying the WKB method. Note that this solution is 
well-known in quantum mechanics as an 'Airy Packet', that spreads without changing its form j53| . 

For ri{x) = kx"^ , the path integral (|4| is quadratic and can be calculated exactly [SO]. The solution is given by 

/K\i/4 1 (^x'^cos{2{DKf/'^t)\ 

"^""'^^'V^dJ y2^sin(2(i?«)V2ly "''P r 2^ s\n{2{Dny/H) ) ' 



The total population corresponding to this solution, N = 1/ \/ cos[2 (Dk)^/^ t], blows up at the finite time 

= (i?«)-V2 I . 

• For r\{x) = /tsign(a;)|a;|", with < a < 2, the WKB analysis can be carried out again and one finds that the 
total population increases as 

«/..„^2c^a.2.2+a^^^ 2a2r(3/2 + 1/a) 



AT ^ exp f ^ (A{af''D'-Kh^+^) ^ V A{a 



(2-|-a)r(3/2)r(l/a) 



For a — 1/2, we arrive at the exp(t^/^) growth that was obtained for a typical Brownian landscape. For a = 1, 
we recover the exp (re^Dt'^/S) growth of the linear landscape. For < a < 1, we obtain the same behavior as for a 
self-affine disordered landscape with Hurst exponent H = a. 



C. Population Growth in High Dimensions 



When the dimension of the underlying substrate exceeds one, various types of random landscapes can arise and 
there is no single landscape which is as natural as the Brownian landscape in Id. Perhaps the closest analog of the 
one-dimensional Brownian landscape is a random landscape that arises by taking the 'height' variable ?7(R) to be a 
Gaussian free field. This class of random manifolds has been widely studied (see e.g. [54] and references therein). 
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The Gaussian free field on a discrete lattice of dimension d is defined by assigning the Gaussian probability to the 
configuration ?y(R): 

exp {- 5][7?(R ± e„) - ry(R)]2} , (62) 

where the sum runs over R e Z"* and over Gq, which arc the unit vectors in the d possible directions. In the continuum 
limit, the summation is replaced by integration, / (Vry)^ dR. Further, one has to introduce a small-scale cut-off to 
regularize correlation functions on short distances [M]. The Gaussian free field is characterized by the following 
height-correlation function: 

{[r]{R) - ri{0)f) - Ini? ford ==2, (63) 
([?7(R) - 77(0)]2) - finite ford > 2. (64) 

One can also define random landscapes by using the Edwards- Wilkinson growth process [3B] which has the same 
2-point correlation functions as above. 



1. Typical Growth 

We estimate the typical growth using the same argument as in one dimension, namely by balancing the typical 
maximal value of the noise over a radius R with the probability for the particles to visit a droplet of size R: 



ntyp exp 



tVlnM- 



R^ 



(65) 



Maximizing with respect to R gives R - t{lnty/^ leading to 



Hyp 



(66) 



In dimensions strictly higher than 2, if we consider a random landscape generated by a Gaussian free field then the 
landscape is statistically fiat, i.e., its width does not vary with the macroscopic length scale R (it depends in fact on 
the microscopic cut-off). Therefore, the typical population grows exponentially. 

To summarize, we have the following growth laws for the typical population under a random landscape 



ntyp 



e*^^^ for c? = 1 



for d = 2 
for d > 3 . 



(67) 



Finally, we note that random surfaces characterized by the logarithmically growing height-correlation function ( 63 ) 
are almost flat. One can consider random surfaces with algebraic height-correlation function (see e.g. (55j ) given by 

([ry(R) — ?7(0)]^) ^ i?''. We then obtain rityp ~ exp 



t i?C/2 _ 



from which we conclude that 



hyp 



exp t 



f(4+C)/(4-C) 



(68) 



2. Average Growth 

The typical growth laws obtained above correspond to a given realization of the Gaussian free field »7(R). The 
average growth of the population is again dominated by the very rare fluctuations of the landscape. To derive the 
average growth we need, as in one dimension, to keep two free parameters: the size R of the droplet visited by the 
particles and the maximum value M reached by ?y(R) over this droplet. 

In two dimensions, a heuristic estimate is found by taking a Gaussian distribution for M of variance Ini?: 



(n) ~ exp 



tM ■ 



2 Ini? 



R^- 

At 



(69) 
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Growth in various d 


Typical 


Average 


One dimension 






Two dimensions 






Higher dimensions {d > 2) 







TABLE II: Asymptotic growth laws in a Gaussian random landscape. In order to regularize short scale fluctuations in d > 2, 
the Gaussian free field is defined on a lattice. For d — 1, the results were already given in Table |l] 



Maximizing with respect to M and i?, we find 

(n)-e''i"*. (70) 

A rigorous derivation of this estimate is perhaps a challenging problem. Nevertheless, we can justify why the distri- 
bution of the maximum value M of a Gaussian free field has a Gaussian tail. Let us first revisit the d = 1 case. The 
full distribution of the maximum of a Brownian process is of course well-known |49j , yet we want to deduce its tail 
in a simple way that will admit a generalization to the Gaussian free field. The tail of the maximum distribution can 
be retrieved by the following simple reasoning: let us consider a Brownian path h{x) of length L with ft,(0) = and 
let M be its maximum value; when M is large, the maximum must be reached in the vicinity of the end of the path 
and such a path contributes by a weight of exp(— | J {h' {x))"^ dx) . The optimal path that has the largest weight is 

therefore obtained by minimazing the integral (h' (x))^ dx . This gives h" — 0, which in conjunction with h{0) — 
and h{L) = M, leads to h{x) = Mx/L. Substituting this expression in the exponential weight, we arrive at the tail 
exp(-M^/2L) which we used in Eq. ( [40| . 

Turning to two dimensions, let us consider a Gaussian random surface h{x, y) over a circular disk of radius R. As 
in d = 1 we expect the maximum value M of the height h to be reached on the rim of the disk. Supposing that the 
optimal surface is rotationally invariant, its statistical weight is given by exp(— ^ Ji^)^'fdr) in radial coordinate r. 
Optimizing this weight with the constraints h{0) = and h{R) = M leads to the Euler-Lagrange equation: 

dh^, 
r—\=0. 



dr V dr 

The solution to this equation is h{r) = and the weight of this optimal path is given by 

/ M2 dr\ ( AP 



exp m9 / — = cxp 



2(lni?)2 r I ' \ 21ni? 

where a short-length scale cut-off (we set it to unity) allows to avoid the small r divergence. This justifies the 
expression of the distribution of M used in equation ( [69| ). 

A similar reasoning can be carried out in higher dimensions. In order to estimate the average growth we again 
need to know the distribution of the maximum. Proceeding as above we find that in three dimensions, the spherically 
symmetric optimal landscape satisfies (''^^) — with respect to the radial coordinate r. This leads to h{r) = 
~^ M(1 — r~^) (taking again the microscopic cut-off to be unity). The corresponding weight behaves as 

exp(— M^/2); there is no dependence on R because the interface is flat at large scales. Using this expression for the 
maximum distribution, we obtain the average growth of the population 



M2 i?2 

tM 

2 At 



[n) ^ exp 

in three dimensions. Maximizing with respect to Al we arrive at 

(n) ^ e*^ . (71) 
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Note that we do not need to optimize with respect to R because the distribution of M, for large values of M, is 
independent of R. A similar reasoning can be carried out in higher dimensions and we find that the behavior is also 
given byjTlJ for all d > 2. 

Table [Il[summarizes our results for the typical and the average population growth in the situation when the random 
landscape is described by a Gaussian free field. 



IV. CONCLUSION 



In this work, we have investigated the evolution of a population of non-interacting particles that undergo diffusion 
and birth/death. The latter depends on the environment, for example, the distribution of nutrients that defines a 
landscape which is assumed to be stationary. Different statistical properties of this landscape lead to a number of 
laws for the growth of the population. In most of the cases, the total population increases in a faster-than-exponential 
manner. This behavior is due to two features of the noise, namely its multiplicative nature and the lack of upper 
bound. Another striking feature is the huge difference between typical and average behaviors. In order to determine 
the average growth law, one has to consider all possible realizations of the random landscape and the average is 
dominated by very rare configurations. Thus extremal statistics play an important role. Some of our analysis relies 
on heuristic arguments. In one dimension we have performed asymptotically exact calculations in the situation where 
the noise is described by a Brownian process. We have shown that the determination of the average population growth 
reduces to calculating a first passage exponential functional of the Brownian motion. This problem can be solved by 
using the Backward Fokker-Planck equation. The asymptotically exact results agree with heuristic predictions. 

Although the basic stochastic differential equation ([I]) has been studied for almost thirty years, there are still many 
open problems. Even in one dimension, it would be interesting to generalize the quantitative approach applicable to 
the Brownian landscape to other random landscapes such as those generated by a fractional Brownian motion. In 
higher dimensions {d > 2), little is rigorously and/or exactly known. To appreciate the challenge, one can think of 
the somewhat related problem of the localization of a quantum particle. 

On a more practical side, it could be interesting to study the transient regime when the system is still far from the 
final asymptotic stage. In the long time limit, non-linear saturation effects, that are ignored in our model, can start 
playing a prominent role. Adding nonlinear terms to the basic equation ([T]) will totally modify its properties in the 
asymptotic regime. This is a challenging mathematical problem that deserves further analysis. 

The work of PLK has been supported by NSF Grant No. CCF-0829541. We are thankful to M. Bauer, F. David, 
B. Duplantier, S. Mallick, and S. Redner for suggestions, help, and encouragement. 
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